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Motivated by the recent observations of incommensurate magnetic order and electric polarization 
in YBaCuFeOs up to temperatures Tn 2 as high as 230K— we report here for the first time a 
model for the incommensurate magnetic structure of this material that we complement with ab- 
initio calculations of the magnetic exchange parameters. Using neutron powder diffraction we show 
that the appearance of polarization below Tn 2 is accompanied by the replacement of the high 
temperature collinear magnetic order by a circular inclined spiral with propagation vector ki=(l/2, 

1/2, l/2zbg). Moreover, we find that the polarization approximately scales with the modulus of 
the magnetic modulation vector q down to the lowest temperature investigated (3 ~K). Further, 
we observe occupational Fe/Cu disorder in the FeOs-CuOs bipyramids, although a preferential 
occupation of such units by Fe-Cu pairs is supported by the observed magnetic order and by density 
functional calculations. We calculate exchange coupling constants for different Fe/Cu distributions 
and show that, for those containing Fe-Cu dimers, the resulting magnetic order is compatible with 
the experimentally observed collinear magnetic structure (kc=(l/2, 1/2, 1/2), Tn 2 > T > Tati = 
440K). Based on these results we discuss possible origins for the incommensurate modulation and 
its coupling with ferroelectricity. 

PACS numbers: 71.27.+a, 71.30.+h, 71.45.Lr, 61.05.fm 


I. INTRODUCTION 

The discovery of materials with strongly coupled mag¬ 
netic and ferroelectric orders has raised a great deal of 
interest in view of their possible use in magnetoelectic 
device applications. In many such materials the sponta¬ 
neous appearance of electric polarization (P) is linked to 
the onset of incommensurate (ICM) magnetic order 
This is often the signature of competing magnetic in¬ 
teractions — and hence is characterized by low ordering 
temperatures. As a result, their promising technological 
multifunctionalities such as the control of magnetism by 
applied electric fields usually occur at temperatures too 
low for most practical applications. 

To date there are only two known examples of switch- 
able, magnetism-driven ferroelectricity at zero field and 
temperatures above 200K. One is cupric oxide CuO^ii^ 
where the low monoclinic symmetry promotes frustrated 
magnetic interactions. A spiral magnetic multiferroic 
phase results in the limited temperature range of 213 
to 230 Kj^ The second example, which has received much 
less attention due to its unavailability as a single crystal 
or thin film, is the layered perovskite YBaCuFeOs. This 
material was first synthesized in 1988^ one year after 
the discovery of high-temperature superconductivity in 
YBa2Cu306+cc^^ Although it is not superconducting^ 
YBaCuFeOs displays magnetism-driven ferroelectricity 


at unexpectedly high temperatures (T < Tn 2 ^ 230K)i 
and in a temperature range more than 10 times larger 
than CuO. The reported P values are also fairly large, 
reaching 0.4 fiC/cw? in powder samples^ In contrast to 
CuO, the ICM magnetism observed in the ferroelectric 
phase of YBCFO is somewhat surprising as the high sym¬ 
metry tetragonal perovskite structure of this material is 
a priori not prone to magnetic frustration. A further 
unanswered question is the nature of the ICM magnetic 
order, which to our best knowledge has not been reported 
although its existence has been known since 1995 To 
understand its origin and the consequent multiferroism, 
a detailed knowledge of the magnetic structure and ex¬ 
change interactions is clearly required. 

Here we report novel neutron diffraction results that 
enable us for the first time to propose a model for the 
ICM magnetic structure of YBaCuFeOs. We also present 
experimental evidence suggesting the existence of a par¬ 
ticular kind of Cu/Fe chemical disorder characterized by 
the existence of Cu-Fe bipyramidal dimers. Density func¬ 
tional theory (DFT) calculations carried out for vari¬ 
ous Fe/Cu distributions provide additional support for 
this scenario. We also calculate the magnetic exchange 
coupling parameters between next and (selected) next- 
nearest neighbours in order to gain insight into the origin 
of the magnetic order and the magnetic frustration in this 
material. To conclude, we discuss the polarization direc- 
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tion and the possible origins of the multiferroicity based 
in the symmetry of the ICM magnetic structure and in 
the similar temperature dependences of the polarization 
and the magnetic modulation parameters. 


II. EXPERIMENTAL DETAILS 


The YBaCuFeOs ceramic sample used in this work was 
prepared by solid state synthesis. Stochiometric amounts 
of Y2O3 (previously pre-annealed at 950°C for 10 hours), 
BaCOs, CuO and Fe2 03 were thoroughly mixed and fired 
at 1050°C in air for 50h. The resulting black powder 
was ground and pressed into pellets which were annealed 
again under the same conditions. The phase purity was 
checked by laboratory x-ray powder diffraction (Brucker 
D8 Advance, Cu Kq,), which indicated the absence of 
foreign phases as well as an excellent crystallinity. The 
cationic distribution in the sample, as determined by x- 
ray fluorescence spectroscopy, was found to be homoge¬ 
neous within a 30//m scale with Y/Ba/Fe/Cu ratios in 
excellent agreement with the nominal composition. The 
oxygen content, as determined from thermogravimetric 
analysis, was 4.95(2). 

DC magnetic susceptibility measurements were carried 
out on a commercial Physical Properties Measuring Sys¬ 
tem (Quantum Design). A cylinder-shaped YBaCuFeOs 
pellet (D = 5mm, H = 7mm) from the same batch as 
the sample used for the neutron diffraction measurements 
was cooled in zero field down to 1.8K. The magnetization 
of the sample was then measured in a magnetic field of 
IT while heating at a constant rate of 2K/min. 

The electric polarization was determined from pyrocur- 
rent measurements which were carried out using a Kei- 
thely 6517B electrometer. A thin pellet (D = 11mm, 
H = 1mm) was sputtered with gold on both faces and 
mounted on the stick of a He cryostat. The sample was 
cooled from RT down to 3K with an electric field of ± 
300V applied between the gold-covered faces. At 3K the 
field was removed and the stray charges (if any) recorded 
during 600s. The pyrocurrent was then measured by 
heating the sample at 20K/min. 

Neutron powder diffraction (NPD) measurements were 
carried out at the Swiss Neutron Source SINQ of the 
Paul Scherrer Institut in Villigen. Several patterns were 
recorded between 1.5 and 500K at the powder diffrac¬ 
tometers HRPT (A=1.1546A)i^ and DMC (A=2.45A).i^ 
The two series of experiments were carried out consec¬ 
utively using the same cryofurnace, whose contribution 
to the background was minimized using oscillating radial 
collimators. All data were analyzed using the Rietveld 
package FullProf SuiteJ^ 


III. CRYSTAL STRUCTURE AND FE/CU 
DISTRIBUTION 

The crystal structure of YBaCuFeOs is displayed in 
Fig.m The tetragonal unit cell can be described as an or¬ 
dered array of layers containing the large Ba^+ ions plus 
two corner-sharing square pyramids separated by Y^+ 
sheets. Equal amounts of Fe^+ and Cu^+ sit inside the 
pyramids, though at different distances from the basal 
plane If the two ions are equally distributed among the 
pyramids as shown in Figs. [1^ and b, the average struc¬ 
ture is centrosymmetric with space group (SG)P^/mmm. 
If the Fe/Cu distribution is asymmetric as in Figs. [It and 
d, the mirror plane containing the Y^+ ions is lost and 
the structure is non-centrosymmetric (SG P^mm). Per¬ 
fect Fe/Gu order along the c axis is a particular case of 
this scenario (see Fig. m). 

Due to the relatively small difference between the Fe^+ 
and Gu^+ ionic radii (0.07Ai^), the occupation of the 
two sites is in fact strongly dependent on the preparation 
method and the two space groups have been reported by 
different authors in the past j^^i^^i^^~ — This has an im¬ 
pact on the magnetic interactions and on the value of 
Tn 2 ^ which displays an important dispersion in the lit¬ 
erature (180 to 24QK)ji^ d^d^i^^~ — A proper comprehen¬ 
sion of the magnetism in YBaGuFeOs thus requires thus 
information about the Gu/Fe distribution within the 2 
available square-pyramidal sites. 

The results obtained for the refinement of the high 
resolution NPD data measured on HRPT (A = 1.1546 
A) at room temperature (RT) are summarized in Tab.lH 
We used two different Fe/Gu distribution models for each 
SG. For P4/mmm we compared a single site occupied by 
50% Fe and 50% Gu with a model with split Fe and Gu 
sites, each of them also with half occupation. As shown 
in Table [H the fit obtained using split sites results in 
significantly better reliability factors. 

A similar conclusion was derived from the fits using 
which, in contrast to centrosymmetric P^/mmm, 
allows us to refine the Fe/Gu occupation. For this space 
group we compared a model with full Fe/Gu order with 
another with split Fe/Gu sites. The agreement between 
the observed and the calculated patterns was again sig¬ 
nificantly better for the second case, which displayed the 
best reliability factors of the four models presented in Ta¬ 
ble |I| (see also Fig. [2]). The occupation of the split Fe/Gu 
sites obtained with this last model is slightly asymmetric 
with approximately 1/3 Fe 2/3 Gu {z ^ 0.25) and 2/3 
Fe -h 1/3 Gu {z ^ 0.75). This indicates that the material 
is non-centrosymmetric on average with a polar axis par¬ 
allel to the 4-fold axis along c. As we will show later, no 
permanent polarization is observed above Tn 2 = 200K, 
probably because the Fe/Gu disorder prevents the exis¬ 
tence of coherence between the Fe/Gu displacements. 

The existence of disorder is also supported by the 
Fourier difference maps calculated for the four models, 
which are displayed in Fig.lH We used the FullProf Suite 
program Fourier— and structure factors corresponding 
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MODEL 1 MODEL 2 MODEL 3 MODEL 4 

T = 300K P4/Tnmm P4/Tnrara P4Tnra P^Tnm 

(1 single Fe/Cu site) (Fe/Cu sites split) (Fe/Cu fully ordered) (Fe/Cu split, partial order) 


a(A) 

3.87323(1) 

3.87326(1) 

3.87325(1) 

3.87325(1) 

c(A) 

7.6651(3) 

7.6655(3) 

7.6651(3) 

7.6655(3) 

Ba 

la (0 0 0) 

la (0 0 0) 

la (0 0 z) 

la (0 0 z) 

z 



0 

0 

t/ll(A2) 

0.0042(4) 

0.0043(3) 

0.0024(4) 

0.0035(3) 

U33{k^) 

0.0189(12) 

0.0231(12) 

0.014(2) 

0.024(2) 

Y 

lb (0 0 1/2) 

lb (0 0 1/2) 

la (0 0 z) 

la (0 0 z) 

z 



0.4931(10) 

0.5053(16) 

Uuo{k^) 

0.00263(17) 

0.00250(14) 

0.00281(19) 

0.00249(18) 

Cu 

2h (1/2 1/2 z) 

2h (1/2 1/2 z) 

lb (1/2 1/2 z) 

lb (1/2 1/2 z) 

Oee 

0.5 

0.5 

1 

0.703(2) 

z 

0.26729(13) 

0.2833(3) 

0.7155(10) 

0.2856(7) 


0.004355i) 

0.0015652) 

0.00402(12) 

0 . 0014953 ) 

Fe 

2h (1/2 1/2 z) 

2h (1/2 1/2 z) 

lb (1/2 1/2 z) 

lb (1/2 1/2 z) 

Oee 

0.5 

0.5 

1 

0.297(2) 

z 

0.26729(13) 

0.2544(2) 

0.2511(10) 

0.2516(8) 


0.004355i) 

0.0015652) 

0.00402(12) 

0.00149(13) 

oi 

Ic (1/2 1/2 0) 

Ic (1/2 1/2 0) 

lb (1/2 1/2 z) 

lb (1/2 1/2 z) 

Z 



-0.0096(14) 

0.0179(15) 

t/ii(A2) 

0.0074(4) 

0.0067(4) 

0.00845) 

0.0067(4) 

t/33(A2) 

0.0136(9) 

0.01455) 

0.0157(11) 

0.0090(17) 

02 

4i (0 1/2 z) 

4i (0 1/2 z) 

2c (1/2 0 z) 

2c (1/2 0 z) 

Z 

0.31603(12) 

0.31601(10) 

0.3077(12) 

0.3265(15) 

f7ii(A2) 

0.0050(3) 

0.0043(3) 

0.0052(3) 

0.0043(3) 

C/22(A2) 

0.0022(3) 

0.00245) 

0.00245) 

0.00235) 

U33{k^) 

0.0113(4) 

0.0090(3) 

0.0112(10) 

0.0086(4) 

02’ 



2c (1/2 0 z) 

2c (1/2 0 z) 

z 



0.6758(12) 

0.6947(15) 

c/ii(A2) 



0.0052(3) 

0.0043(3) 

f722(A2) 



0.00245) 

0.00235) 

t/33(A2) 



0.0112(10)) 

0.00865) 

Chi^ 

2.64 

2.01 

2.46 

1.98 

Rp 

4.52 

3.90 

4.33 

3.86 

Rwp 

5.74 

5.01 

5.54 

4.96 

Ruragg 

5.34 

4.05 

4.98 

3.81 


TABLE I: Atomic coordinates, thermal parameters and Fe/Cu mag¬ 
netic moments of YBaCuFeOs at 300K, as refined in the space groups 
P4/mmm and Pfmm (both with Z=l) using the neutron powder diffrac¬ 
tion data recorded on HRPT (A = 1.1546 A). The reliability factors of 
the different models are also provided. 


to reciprocal space vectors H with modulus smaller Fourier trausforms of the differeuce betweeu the observed 

thau 10.7 A“^. The coutour plots iu Fig. [T]represeut the aud the calculated ueutrou scatteriug deusity at RT. The 
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Ba^ 


a) P4/mmm (Model 1) 

— ® 



L 
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b) P4/mmm (Model 2) 

Fe 



-0.5 


Cu 


0.5 


Cu 

Fe, 
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0.5 




d) P4mm (Model 4) 



FIG. 1. (Color online) Models for the crystal structure of YBaCuFeOs used for the Rietveld fits of the RT NPD data recorded 
on HRPT (see also Table |T|. (a) Centrosymmetric (P4/mmm) with a single Fe/Cu site, (b) Centrosymmetric (P4/mmm) with 
Fe/Cu sites split, (c) Non-centrosymmetric (P4mm) with perfect Fe/Cu order along c. (d) Non-centrosymmetric (P4mm) with 
partial Fe/Cu order. The contour plots are Fourier difference maps of the (x, 1/2, z) plane showing the neutron scattering 
density not reproduced by each of the models. 


results obtained for the two ordered models (1 and 3 in 
Table [I]) show the existence of scattering density (bright 
yellow spots) not reproduced by these models above and 
below the refined Cu/Fe positions (black circles). These 
spots are absent for the disordered models 2 and 4, with 
results slightly better for the second one. The models in¬ 
volving Cu/Fe disorder (that we assume to be random in 
our fits) provide thus a better description of the experi¬ 
mental data, even if the presence of correlations between 
the Fe and Cu site occupations cannot be completely dis¬ 
regarded. In particular, it has been suggested that the 
bipyramidal units linked by the apical oxygen 01 could 
always host Cu-Fe pairs, which would be randomly ar¬ 
ranged in the structureSince the Fe-01 and Cu-01 
apical distances are slightly different, such a distribution 
is expected to give rise to lower microstrains along c than 
one with coexisting Cu-Cu (long), Fe-Fe (short) and Cu- 
Fe (intermediate) pairs. We show later that both the 
proposed magnetic structures and the DFT calculations 
favor such a scenario. 


IV. MAGNETIC TRANSITIONS 

Fig. [3^ shows the temperature dependence of the low- 
angle region of the NPD patterns recorded on DMC. Two 
phase transitions involving the appearance of new Bragg 
reflections are clearly observable at T/vi ^ 440K and T/v 2 
^ 200K. The new peaks appearing below Tati correspond 
to the propagation vector kc = (1/2, 1/2, 1/2), indicating 
that the magnetic structure is commensurate (CM) with 
the crystallographic unit cell between 440 and 200K in 
agreement with previous reports 

Below T/V 2 two satellites appear around each CM mag¬ 
netic reflection. They can be indexed with the propaga¬ 
tion vector k^ = (l/2, 1/2, l/2±g), which involves an ICM 
modulation of the magnetic moments along the c axis. 
The modulation parameter q increases continuously with 
decreasing temperature and remains ICM down to 1.5K, 
see Fig.m The onset of the ICM magnetic order coincides 
with a sharp anomaly in the magnetic susceptibility and 
with the appearance of electric polarization P (Fig. 0]), 
indicating a direct relationship between the two phenom¬ 
ena. Below 50K P reaches ^ 0.64/iC/cm^, close to the 
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2Theta (degrees) 


FIG. 2. (Color online) NPD pattern measured on HRPT at 
RT. Red crosses: observed data. Black lines: Rietveld fit 
obtained using the model 4 of Table H] The vertical ticks 
indicate the positions of the Bragg reflections for the crystal 
(upper row) and the commensurate AFM magnetic structure 
with kc = (1/2, 1/2, 1/2) (lower row, see also Fig. [5^). 

value reported by Kundysi (^ 0.4//C/cm^) and about 
10 times larger than the one reported by Kawamura?^ 
(0.04/iC/cm^). The temperature dependence of P closely 
follows that of g, a behavior that we will address later in 
the text. 

Fig.[3]3 shows the evolution of the integrated intensities 
of the first magnetic peak and its satellites. Below Tn 2 
the main CM reflection (1/2, 1/2, 1/2) starts to decrease. 
At the same time the intensities of the two satellites start 
growing. At 1.5K the CM peak is still visible, but its 
intensity is much lower than that of the ICM satellites. A 
similar behaviour is observed for all CM/ICM reflection 
sets, see Figs. [3^ and [6l Such behavior contrasts with 
previous reports, where the ICM satellites were either 
absent^ or much less intense than the CM reflections 
at all temperaturesThis is illustrated in Fig. [3t, 
where the region around the first CM reflection (1/2, 
1/2, 1/2) measured at 1.5K for our sample (black) and 
a sample prepared according to the method described in 
(blue) are displayed. The improved quality of the 
ICM magnetic reflection set obtained with our synthesis 
procedure, together with the larger number of reflections 
measured compared with previous worksenables 
us for the first time to propose a model for the ICM 
magnetic structure of YBaCuFeOs. 


V. MAGNETIC STRUCTURES 

We used representation analysis to find the possible 
magnetic moment arrangements compatible with the the 
space groups P4/mmm and P4mm. The characters of 
the irreducible representations {irreps) of the little group 
Gj^, formed by the operations of the space group G which 
leave the /g-vector invariant, are shown in the Appendix 


a) DMC A = 2.45 A 



T(K) 2Theta (degrees) 


FIG. 3. (Golor online) a) 2D contour plot showing the temper¬ 
ature dependence of the NPD patterns for YBaGuFeOs- The 
patterns at 1.5K and 300K are shown separately, b) Temper¬ 
ature dependence of the integrated intensity of the (1/2, 1/2, 
1/2) magnetic reflection and its incommensurate satellites, c) 
Portion of the NPD patterns showing these reflections for our 
sample (230K and 1.5K) and a sample prepared according to 
ref.— (1.5K). d) Temperature dependence the Fe^^ magnetic 
moment. Inset: angle 0 between the magnetic moment direc¬ 
tion (Tni<T<Tn 2 ) and the spiral plane (T<Tn 2 ) with the c 
axis. 

A and B for both P4/mmm and P4mm and the two prop¬ 
agation vectors kc and k^ together with tables with the 
symmetry-allowed basis functions associated with each 
irrep. Such tables indicate that, if only one irrep be¬ 
comes active below each of the magnetic transitions, the 
direction of the moments is restricted to either the ab 
plane or along the c axis. We find instead that, at all 
temperatures, the best fits to the data correspond to in¬ 
clined arrangements needing the combination of 2 irreps^ 
see the Appendix A and B for details. 

Two collinear models and a commensurate helix are 
compatible with our observations and they give rise to 
the same neutron powder diffraction patterns if Cu and 
Fe are located exactly at 2 ; = 1/4 and 2 ; = 3/4. For 
the refined Fe and Cu positions, which have z coordi¬ 
nates rather close, but not identical to these values (see 
Table [I]) the three models give rise to slightly different 
intensities. As shown in the Fig. [TOl the best agreement 
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FIG. 4. Red dots: DC magnetic susceptibility of YBaCuFeOs 
measured at IT by heating after cooling in zero field. Black 
open squares: Incommensurate modulation vector q (in re¬ 
ciprocal lattice c* units). Blue dotted/continuous lines: nor¬ 
malized electric polarization measured by applying an electric 
field of di 300V. In order to use the same axis as q the polar¬ 
ization values have been normalized to their saturation value 
0.64/iC/cm^) and further divided by a constant value of 
- 10.2 . 


a) 230 K 


b) 1.5 K 




k, = (72 72 72) 


ki = (72 72 72 ± q) 


FIG. 5. (Color online)Magnetic structures of YBaCuFeOs. a) 
Colinear magnetic order at 230K. b) Two views of the circular 
spiral order at 1.5 K. For clarity, only one crystal cell within 
the ab plane is shown. 


corresponds to the collinear magnetic structure displayed 
in Fig. [5^. Spins in the ab plane are antiferromagneti- 
cally (AFM) coupled whereas the alignment along c al¬ 
ternates: it is AFM across the 0-free Y planes and fer¬ 
romagnetic (FM) inside of the bipyramidal blocks. Note 
that the orientation of the spins the ab plane is arbitrary 
because their direction in the plane perpendicular to the 
4-fold axis cannot be determined from NPD. The angle 
0 between magnetic moments and c axis decreases con¬ 
tinuously from 75° (at 400K) to 65° (at 200K), see inset 
of Fig. [^. 

Below Tn 2 two propagation vectors kc and are 
present down to the lowest temperature investigated, see 
Fig- Eb and Fig. [6l The magnetic intensities can be de¬ 
scribed as arising from either a multi-A j arrangement or 
from two distinct magnetic phases; our NPD data alone 
cannot distinguish between the two possibilities. How¬ 
ever, the fact that the intensities of both reflection sets 
have opposite temperature dependencies below Tn 2 fa¬ 
vors the second one. 

As in the CM case, two different models give rise to 
very similar magnetic intensities. However, the excellent 
statistics of the DMC data and the fact that Fe and Cu z 
coordinates differ significantly from z = 0.25 and 0.75 at 
all temperatures allow us to distinguish them, see Fig.[6l 
Fig. [TT] and further details in the Appendix B. The best 
agreement with the experimental data corresponds to the 
circular helix displayed in Fig. [5b - Such a spin arrange¬ 


ment involves the preservation of the AFM coupling be¬ 
tween the 3d metal sites without connecting oxygen ob¬ 
served in the commensurate phase, see Fig. jSb- It also 
implies the loss of the FM coupling within the bipyra¬ 
mids, suggesting that this magnetic coupling is more af¬ 
fected by the thermal evolution of the structure than is 
the one across the 0-free Y layers. 

Our results indicate that the plane of the helix forms 
an angle 0 with the c axis which decreases continuously 
with temperature {0 ^ 65° at 200K, 0 ^ 45° at 1.5K, the 
anomalously low value at 150K is probably due to the 
strong superposition of the 2 incommensurate satellites 
at this temperature). The ICM magnetic order evolves 
thus from an inclined helix towards a cycloid with de¬ 
creasing temperature. The reasons for this behaviour 
are unclear, but it could be related to different Cu^+ 
and Fe^+ magnetic anisotropies and their effect on the 
temperature dependence of the magnetic moment orien¬ 
tations. 

Using the previously described models for the CM and 
ICM magnetic structures we determined their fractions 
in the sample below Tn 2 ‘ By restricting the modulus and 
the 0 angle of the Cu^+ and Fe^+ magnetic moments in 
both phases to be identical and their ratio to be the same 
as for their spin-only values (1 to 5) we obtain 8% (CM) 
and 92 % (ICM). This contrasts with previous studies, 
where most of the sample displayed CM magnetic order 
below T/V 2 and only a small fraction was able to undergo 
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the CM ^ ICM phase transition. As we discuss in the 
next sections, differences in the Fe/Cu distribution are 
the most probable origin of these differences. 

The evolution of the Fe^+ magnetic moment is shown 
in Fig. [3li. Its value at 1.5 K is 3.74(2) (0.748(4)/iB 
for Cu^+), about 1/3 reduced with respect to those ex¬ 
pected for the free-ion, spin-only moments. Such reduc¬ 
tion may, among other reasons, also be related to Fe/Cu 
disorder. The dip around Tn 2 is not an artifact from the 
fits since the same dip is observed in the total (central 
peak -f satellites) integrated intensity as a function of 
temperature (Caignaert and co-workers found a similar 
behaviour, see^^). This suggests that the incommensu¬ 
rate domains are large enough to produce a mesurable 
contribution to satellites only a few degrees below Tn 2 - 



20 40 60 80 


2Theta (degrees) 

FIG. 6. (Color online) NPD pattern measured on DMC at 
1.5K. Red crosses: observed data. Black lines: Calculated 
intensities. The vertical ticks indicate the positions of the 
Bragg reflections for the crystal (upper row), the CM mag¬ 
netic structure with kc = (1/2, 1/2, 1/2) (intermediate row) 
and the ICM magnetic structure with kc = (1/2, 1/2, l/2q=g) 
(lower row). 


VI. MAGNETIC COUPLINGS FROM 
GOODENOUGH-KANAMORI-ANDERSON 
RULES 

To understand the origin of the magnetic structures 
we first consider the Goodenough-Kanamori-Anderson 
(GKA) rules of superexchange. For Fe^+ (3d^) and 
Cu^+ (3d^) in the square-pyramidal arrangement of the 
YBaCuFeOs structure, these rules predict strong AFM 
nearest-neighbour (NN) couplings within the ab plane 
in agreement with our observations. Interestingly, the 
sign of the eouplings is independent of the Cu/Fe distri¬ 
bution. Along the c axis the situation is more complex. 
On one side, superexchange between NN is not possible 
across the Y layers due to the absence of the apical O. 
On the other, the sign of the coupling within the bipyra¬ 
mids is expected to be strongly dependent on the Fe/Cu 
distribution (AFM for two Fe ions, FM for a Fe-Cu pair 


and negligibly small for two Cu ions). This agrees with 
the observation of an ICM modulation only along this 
direction and suggests that the appearance of magnetic 
frustration below TiV 2 could be due to a temperature- 
driven imbalance of the magnetic couplings along the c 
axis. Such a scenario is consistent with the continuous 
evolution of the modulation parameter q with decreasing 
temperature (Fig. S]). Note also that, within the GKA 
framework and for the CM phase, the experimentally ob¬ 
served FM coupling within the bipyramids is only possi¬ 
ble if they are oeeupied by a Cu-Fe pair. 


VII. AB INITIO CALCULATIONS 

Further insight requires detailed quantitative calcula¬ 
tions of the exchange interactions that we obtain using ab 
initio calculations with the Local Spin Density Approx¬ 
imation plus Hubbard U (LSDA-fU)^ to density func¬ 
tional theory as implemented in the Vienna ab-initio sim¬ 
ulation package (VASP) j^^i^^ For these calculations we 
used projector augmented wave potentials with 3d and 3p 
orbitals treated as valence states for both Fe and Cu ions. 
The values of on-site effective Coulomb interaction were 
set to Upe = 5 eV and Ucu = 8 eV while the on-site effec¬ 
tive Hand’s couplings were = 1 eV and Jcu = ^ for Fe 
and Cu, respectively. For all the calculations, except for 
those relative to some of the next-nearest-neighbor mag¬ 
netic couplings along c (see Appendix C), we considered 
a supercell with a = ^/2 die and c = 2cc (ac and Cc are 
the crystallographic unit cell parameters) containing four 
formula units, as depicted in Fig.[3 a F-centered k-point 
grid of size 8x8x4 and an energy cut-off Ecut = 600 eV. 
Relaxations of the reduced ionic positions in the unit cell 
were performed until all atomic forces were below 2 • 10“^ 
eV/A and using the experimental lattice parameters at 
1.5K (a=b=5.462 A, c=15.258 A for the supercell.). 


A. Cu/Fe distribution 

First, we investigated the stability of the crystal struc¬ 
ture of YBaCuFeOs under different Fe^+/Cu^+ distri¬ 
butions. The fractional coordinates of all atoms were 
relaxed starting from different arrangements for the two 
transition metal ions, keeping the commensurate antifer¬ 
romagnetic order of Fig. [5^, fixing the lattice parameters 
to the experimentally observed values at 1.5K and in ab¬ 
sence of the spin-orbit coupling. Fig.^shows the different 
configurations for the Fe/Cu distribution and the energy 
of the relaxed structures compared with the lowest en¬ 
ergy one. We see that there are two energy hierarchies 
depending on which ions are present in the bipyramids. 
Low energy configurations (up to 0.214 eV ^ 2500 K), 
panels (a) to (e) within the blue frame, contain both Fe^+ 
and Cu2+ in all bipyramids, while those containing two 
Cu^+ and/or two Fe^+ ions in at least one of the bipyra¬ 
mids have higher energy (from 0.996 eV ~ 11600 K), 
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panels (f) to (j) within the red frame . 



alts XI H XT H 



0.000 eV 0.015 eV 0.095 eV 0.190 eV 0.214 eV 



nn n nn 



0.996 eV 0.996 eV 0.998 eV 1.558 eV 2.312 eV 


FIG. 7. (Color online) Three-dimensional view of the su¬ 
percell with a = a.c and c = 2cc used in the ab-initio 
calculations (left upper panel) and its projection on the ab 
plane (right upper panel) for one of the Fe^^ /Cu^^ distribu¬ 
tions considered. Note that the a and b axes are rotated by j 
with respect to those in Fig.[T] Panels from (a) to (j) show the 
projection on the ac plane of the considered supercells with 
different orderings of Fe^"^ ions (golden square pyramids) and 
(blue square pyramids) and their energies relative to 
that of structure (a). For clarity Ba^^ and ions are not 
shown. The blue and red frames surround respectively the 
low-energy and high-energy sets mentioned in section VILA. 


B. Nearest-neighbour magnetic couplings and CM 
magnetic order 

Next, we evaluated the exchange coupling constants for 
the obtained relaxed structures of the low energy config¬ 
urations (a) to (e), i.e., those where the bipyramids are 
always occupied by a Fe-Cu pair. We assumed that the 
magnetic interactions are described by the Heisenberg 
Hamiltonian 

H = ( 1 ) 

hj 

where, S-pe = 5/2, 5'cu = 1/2 and label the magnetic 
sites, and performed collinear spin polarized calculations 


in the absence of spin-orbit coupling to extract the ex¬ 
change couplings as described in Appendix C. We note 
that in this minimal model we assume the presence of 
only nearest-neighbour (NN) Heisenberg exchanges and 
do not consider higher order terms such as biquadratic 
and ring exchange which have been found to play a role 
in cuprates22ii2i and, more recently, in manganitesj^ 

Results for the exchange coupling constants are listed 
in Tab.im Remarkably, for any pair of NN, the sign of the 
magnetic interactions is the same for all five configura¬ 
tions. The exchange constant, J_l, between neighboring 
spins in the same tetragonal plane is AFM and is by far 
the strongest coupling for all of configurations. JyQ and 
J||, which correspond to the exchange couplings between 
ions in neighboring planes separated by an oxygen layer 
or not are, respectively, FM and AFM and, generally, 
take values smaller than J±. This confirms the c axis as 
the direction with the weakest couplings, which is the di¬ 
rection where an incommensurate wave vector is observed 
below T 7 V 2 - Interestingly, the only FM exchange con¬ 
stant is J ||05 which corresponds to the coupling within 
the bipyramids occupied by a Cu-Fe pair. 

We note that the signs of all couplings in configura¬ 
tions (a) to (e) are in agreement with the Goodenough- 
Kanamori-Anderson rules as well as with the observed 
CM magnetic structure (Tyri < T < Tn 2 ), see Fig. [5^. 
This is not the case for the higher energy configurations 
(f) to (j), where the sign of JyQ for bipyramids occu¬ 
pied by an Fe-Fe pair is AFM and hence not compatible 
with the observed magnetic orderTo summarize, the 
configurations containing exclusively Cu-Fe pairs in the 
bipyramids are more stable, and the coupling within such 
dimers is the only FM one among NN interactions. Since, 
as shown in Fig. [TOl the magnetic structure which gives 
the best reliability indexes involves FM couplings within 
the bipyramids, the only possibility is that such units are 
occupied by Cu-Fe pairs. This result strongly supports 
the existence of Cu-Fe “dimers” as necessary condition 
to stabilize the observed CM magnetic structure. More¬ 
over, to be consistent with the results displayed in Fig.[Tl 
these dimers should be randomly distributed. 


C. Next-nearest-neighbour magnetic couplings and 
ICM magnetic order 

We will focus now on the ICM magnetic structure ob¬ 
served below T/V 2 - To obtain further insight on its origin 
we considered the effect of next-nearest-neighbor Jnnn 
magnetic couplings along the c direction. The values ob¬ 
tained by ab initio calculations are summarized in the 
fourth column of Tab. HD For each of the configurations 
(a), (b), (d) and (f) in Fig.Othere are only two inequiva¬ 
lent Jatatat’s. For configuration (c) their number is much 
larger (eigth). Therefore for this configuration we ex¬ 
tract only Jsj, Ji^ 5 , J 2,6 and J 2 +c ,6 and assume that 

Js-l-c,? Jl-|-c,5 J4,8 ^J2 -|-c, 6 and J4-|-c,8 ~J2,6* 

To estimate the size of Jnnn necessary to give rise to 
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J± (meV) 

J|l (meV) 

J|lO (meV) 

Jnnn (meV) 

(a) 

Ji 2 = 134.5 
Jr ,8 = 8.7 

Ji ,3 = 10.6 
J 5 ,r = 2.8 

J 5.3 = -1.6 

Ji 5 = -0.05 

Jl + c ,5 = - 0.01 

(b) 

Ji ,2 = 129.9 
J 3,4 = 8.9 

Ji ,3 = 1.4 

J 3.5 = -1.6 

Ji ,5 = 0.07 
J 3 ,r = 0.19 

(c) 

J 3,4 = 28.6 
Ji ^2 = 133.0 
Jr ,8 = 8.7 
J 5,6 = 28.2 

J5,7 = 1.1 

J 2,4 = 8.9 
Ji ,3 = 1.7 
J 6,8 = 3.0 

J 2.8 = -1.5 
Jl,7 = -1.5 
J3,5 = -1.7 
J6,4 = - 1.7 

J 3 ,r = 0.20 
Ji ,5 = 0.09 
J 2,6 = - 0.07 
J 2 +c ,6 = - 0.01 

(d) 

Ji ,2 = 28.3 

Jl ,3 = 3.1 
Jr ,5 = 7.5 

J 3,5 = - 1.6 

Ji 5 = -0.04 
Jl+c,5=-0.01 

(e) 

Ji ,2 = 28.3 

Jl ,3 = 1.3 

J 3.5 = - 1.6 

Ji ,5 = 0.12 
J 2,6 = 0.20 


TABLE IL (Color online) Calculated exchange coupling 
constants for the different configurations in Fig. [T] Posi¬ 
tive/negative signs correspond to AFM/FM interactions, re¬ 
spectively. J± , J||o, J|| and Jnnn are defined in the main 
text. To avoid ambiguities the exchange coupling constants 
Jij are also labeled by the position i,j of the magnetic ions 
as in Fig. [71 Note that according to Eq. [5]the energy scale for 
each coupling is obtained by multiplying it by SiSj, that is: 
^ for Fe^^-Fe^^ coupling (brown), | for Fe^^-Cu^^ coupling 
(black) and j for Cu^^-Cu^^ coupling (blue). 

a magnetic spiral state we calculated the phase diagrams 
for the minimal energy state of Eq. [5] by tuning the in¬ 
equivalent J AT AT AT’s of cach Configuration between -1 and 
1 meV and under the assumption that the magnetic order 
is of spiral type. Calculations for configuration (c) were 
not considered due to the large number of inequivalent 
next-nearest-neighbor couplings along c. We assume the 
magnetic structure to be described by a spiral ansatz of 
the form 


S^(Ri) = cos{(f>^ + q ■ Ri)vi + + q ■ Ri)v 2 , (2) 

where vi and V 2 are two orthogonal unit vectors, /i = 
1,..., 8 labels the magnetic ions in the unit cells of Fig. 0 
q is the wave vector, (j)^ is the angle of the /i-th spin inside 
the unit cell and is the position of the i-th unit cell. 
Inserting Eq. [2] in Eq. 1, considering the exchange of the 
first three columns of Tab.lIIland minimizing numerically 
the obtained energy with respect to q and we obtained 
the states depicted in the phase diagrams Fig. [8] as a 
function of inequivalent Jnnn- 

As shown in Fig [51 incommensurate states are stable 
only if at least one of the NNN couplings is ferromagnetic 
(i.e., if it frustrates the collinear commensurate order). 
According to Tab. im the only ferromagnetic NNN cou¬ 
plings are those between a Cu^+ and a Fe^+ which are 
present in configuration a) and d). However, the calcu¬ 
lated strength of these couplings (green dots) is too small 
to stabilize a spiral state. Therefore, the phase diagrams 
in Fig. [8] obtained with the ansatz Eq. [2] indicate that 
the next-nearest-neighbor couplings obtained in Tab. [Ill 
are either too weak or of the wrong sign to stabilize a 
spiral state. However, due to the uncertainty present in 



;,UmeV| y.^|mcV| 

FIG. 8. (Color online) Phase diagrams obtained using the 
spiral ansatz Eq. [2] for configurations (a), (b), (d) and (e). 
The blue color indicates the commensurate collinear structure 
shown in the left pannel of Fig. 3. The red color and the 
yellow color indicate, respectively, an incommensurate spiral 
state and a commensurate non-collinear state. The green dot 
indicates the values of the NNN in Tab. 1. 

the ab-initio calculations of exchange couplings, we can¬ 
not completely exclude NNN coupling as the origin of the 
incommensurate state. 


VIII. POSSIBLE ORIGIN OF ELECTRIC 
POLARIZATION 

To finish, we briefly address the questions of the direc¬ 
tion of the electric polarization in the incommensurate 
phase and the origin of the magnetoelectric coupling. Al¬ 
though they can not be properly answered due to the lack 
of single crystals, it is possible to make some predictions 
using the symmetry properties of the refined ICM mag¬ 
netic structure and the various mechanisms proposed in 
the literature. 

Looking at Fig. [Sb, we see that the only symmetry 
element of the depicted inclined spiral is a binary axis 
along the b direction. Its point group is thus . 2 ., which 
is compatible with the existence of P only along b. Inter¬ 
estingly, the ionic, non-switchable polarization expected 
in the case of a perfect Cu/Fe order would appear along 
c. 

This prediction agrees with those of the Landau The¬ 
ory formalisms proposed in Ref. [fil as well as in Ref. Isol . 
where the magnetoelectric coupling is described by a tri- 
linear interaction term. According to the formalism in 
Ref.S P transforms as the product of the two irreps 
describing the magnetic structure. Looking at Tables 0 
and IVH[ it is easy to see that P is not allowed for kc 
but it may exist in the ICM phase only within the ab 
plane. 
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Similar conclusions can be derived from the 
Dzyaloshinskii-Moriya (DM) and/or spin-current 
models for magnetoelectric coupling. According to these 
mechanisms, the polarization direction is expected to 
be given either by eij x Si x Sj or g x Si x Sj, 
were Cij is the vector connecting the i and j sites and 
q = (0, 0, kz) is the magnetic modulation vector. In 
the incommensurate phase q || eij, so the two models 
make identical predictions. From Fig. [Sb, we see that in 
both cases the polarization is expected to be along h, in 
agreement with the conclusions derived from symmetry 
arguments. 

We also speculate that the similar temperature depen¬ 
dence of P and q shown in Fig. 0] might indicate that 
Dzyaloshinskii-Moriya (DM) and/or spin-current mech¬ 
anisms could be responsible for the ferroelectricity. This 
may look surprising in view of the large value of the polar¬ 
ization in YBaCuFeOs (0.64 ^ /iC/cm^). It is however 
not unreasonable since TbMnOs, where multiferroicity 
is believed to originate from these mechanisms, displays 
a polarization which is only 6 times smaller ( ^ 0.09 
/iC/cm^). 

From the above mentioned mechanisms, one would ex¬ 
pect P oc c X X Sj, where i and j are nearest-neighbor 
magnetic sites along c. As shown in Fig. [3]1, the size of 
magnetic moments below T/v 2 is approximately constant. 
Therefore, one can assume the magnetic spiral state de¬ 
scribed by Eq. [21 where /i = 1,2 labels the magnetic 
sites in Fig. (H vi = sin(6>)a -h cos(6>)c, V 2 = sin(6>)b and 
the spin value is constant. Moreover, as the magnetic 
moments of neighboring magnetic ions along 2 ; not sepa¬ 
rated by oxygen are antiparallel (see Fig.jSb and Table (III 
), one can assume 02 = 0i +7r. Under these assumptions 
the electric polarization given by DM and/or spin-current 
coupling is P 5 = Ccos(6>) sin(27r(^ — q)) where C is a 
temperature-independent constant. 

Fig. [9] shows the comparison between the measured val¬ 
ues of P (blue dots) and those obtained from the previ¬ 
ous expression setting C in such a way that at the lower 
temperature the calculated polarization coincides with 
the observed value. The points represented by black dots 
were obtained using the measured values of q and 0 while 
for the points represented by red diamonds the values of 
the linear interpolation of 0 below 190 K were used (see 
inset of Fig. Taking into account that the only free 
parameter is the proportionality constant C, the agree¬ 
ment between the observed and calculated temperature 
dependence of P is reasonably satisfactory. 

According to the scenario described above, the ap¬ 
proximate proportionality of P and q originates from 
the relatively small changes in q (which implies oc 
27rCgcos(6>)) and the fact that 0 does not change dra¬ 
matically below Tn 2 - We also note that this behavior 
would not hold if the size of magnetic moment would 
strongly depend on temperature below Tn 2 , e.g., if the 
magnetic spiral phase would appear close to Tati- 



FIG. 9. (Color online) Temperature dependence of the elec¬ 
tric polarization (blue dots) and estimations using 0 and q as 
described in the text (red diamonds and black open dots). 

IX. SUMMARY AND CONCLUSIONS 

As mentioned in the introduction, YBaCuFeOs is 
together with CuO the only known material to dis¬ 
play switchable, magnetism-driven ferroelectricity above 
200K. A prerequisite to understand the existence of mul¬ 
tiferroicity at such high temperatures in YBaCuFeOs is 
the knowledge of its ICM magnetic structure. However, 
it was never reported even if its existence is known since 
1995. In this study we have successfully synthesized 
YBaCuFeOs ceramic samples of unprecedented quality 
and we have conducted new, high-resolution and high- 
intensity neutron diffraction measurements that enabled 
us for the first time to solve the incommensurate mag¬ 
netic structure of YBaCuFeOs. Our results are consistent 
with the replacement of the collinear magnetic order ex¬ 
isting between Tn 2 and Tati = 340K by a circular spiral 
with temperature dependent inclination described by the 
propagation vector ki=(l/2, 1/2, l/2Pq). 

The origin of the ferroelectricity observed below T/v 2 
and its coupling with the incommensurate magnetic or¬ 
der can not be fully explained from our results alone, 
in particular because our ceramic samples do not allow 
to determine the direction of the polarization. However, 
the symmetry of the observed magnetic spiral suggest 
that P should be within the ab plane. Also, the very 
similar temperature dependence of the magnetic modu¬ 
lation parameter q and the polarization indicate that the 
Dzyaloshinskii-Moriya and/or spin-current mechanisms 
could be at the origin of the magnetoelectric coupling. 
Further experimental work, preferably on single crystals, 
will be necessary to get additional insight about the va¬ 
lidity of this scenario. 

We have also investigated the crystal structure, in par¬ 
ticular the Fe/Cu distribution between the two available 
square-pyramidal sites, and we find clear evidence for the 
existence of occupational disorder. The observed CM 
magnetic structure suggests however that the bipyrami- 
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dal units may be preferentially occupied by Cu-Fe FM 
pairs. This finding is supported by ab initio calculations, 
which suggest that Cu/Fe distributions containing ex¬ 
clusively Cu-Fe dimers have lower ground-state energies 
than those where Cu-Fe, Fe-Fe and/or Cu-Cu pairs co¬ 
exist. Moreover, the exchange parameters calculated for 
different Fe/Cu distributions are compatible with the ob¬ 
served CM magnetic order only if the bipyramidal units 
are occupied by Cu-Fe pairs. 

To summarize, all our experimental and theoretical re¬ 
sults converge towards the same scenario, i.e., that the 
existence of Cu-Fe dimers may be a necessary condition 
for the existence of the CM magnetic structure. Elastic 
diffuse scattering, which provides information about oc¬ 
cupational correlations, could be used as a further check 
if single crystals become available in the future. 

As in the case of CuO^ the observation of incom¬ 
mensurate magnetic order below temperatures as high as 
200K is most probably related to the large values of some 
of the exchange constants. However, in our case, it is 
also linked to the presence of a crystallographic direction 
where the magnetic couplings are particularly unstable. 
As seen in Fig. 5a, the CM magnetic order can be seen as 
an array of bilayers containing weakly coupled Cu-Fe FM 
bipyramidal dimers coupled through weak AFM interac¬ 
tions along the c axis. These couplings might be frus¬ 
trated by several mechanisms such as NNN interactions, 
bond disorder or higher order magnetic couplings which 
could eventually give rise to a magnetic spiral. We have 
briefly investigated the first possibility and we find that 
the calculated NNN magnetic exchanges along the c axis 
are either too small or have the wrong sign to stabilize 
a spiral state. Although we do not completely exclude 
this mechanism due to the uncertainty in the ab-initio 
calculations for small exchange couplings, investigation 
of the two remaining options is presently in progress. 

The calculated exchanges highlight the strong depen¬ 
dence of the magnetic interactions on the Cu/Fe distri¬ 
bution. This could explain the huge dispersion in the 
Tjsfi and Tn 2 values reported in previous studies, where 
different preparation synthesis procedures were used. In 
general, methods favouring an intimate mixing of the 4 
involved cations seem to result in higher Tn 2 and larger 
fractions of incommensurate phase below Tn 2 - Although 
we don’t have any experimental proof, we speculate that 
such methods may enable the material to reach more eas¬ 
ily the most energetically favourable Cu/Fe distribution, 
i.e., that with bipyramids populated exclusively by Cu- 
Fe dimers. Because in this case all couplings along the 
c axis are weak, it may be easier to frustrate them and 
to generate incommensurate magnetic order. Systematic 
investigations of this matter are currently in progress. 

As a concluding remark, we would like to point out 
that other perovskite oxides have been reported to dis¬ 
play the same crystal structure as YBaCuFeOs. How¬ 
ever, the incommensurate magnetic order and the con¬ 
comitant spontaneous polarization have been only found 
in two other members of the RBaCuFeOs family (R = 


Tm and Lu)#^ The reasons behind this are not yet fully 
understood, but they probably stem from the particu¬ 
lar combination of exchange constants associated to the 
presence of Fe^+ and Cu^+, especially to the weak, alter¬ 
nating FM and AFM couplings along c axis. Since the 
magnetic order is particularly unstable along this direc¬ 
tion, the application of external perturbations may pro¬ 
vide further insight about the hierarchy of the involved 
magnetic interactions. Tuning such weak couplings along 
c (and hence Tn 2 ) could be achieved either by chemical 
substitution, by the application of external pressure or, 
in view if the layered nature of YBaCuFeOs, through 
the fabrication of heterostructures. Given the high value 
of T/V 2 , this qualifies YBaCuFeOs as one of the most 
promising structural frameworks to search for room tem¬ 
perature multiferroics. 


X. APPENDIX 

The possible magnetic structures compatible with the 
space groups PA/mmm and PAmm were determined 
using representation analysis. The only magnetic ions 
in YBaCuFeOs are Fe^+ and Cu^+, and there is only 
one Fe and one Cu atom per crystallographic unit cell. 
The labeling of the different sites for both space groups 
is displayed in Tab. IHIl The main difference between 
P4/mmm and PAmm is that the sites 1 {z ^ 0.25) and 2 
{z ^ 0.75) are symmetry-related in the first one whereas 
they are not in the second. Since for a given space group 
Cu and Fe occupy identical Wyckoff positions, we discuss 
only the Fe case in the following. All conclusions apply 
also to Cu. 

TABLE III. Atomic positions of Fe and Cu within the unit 
cell in the space groups P4/mmm (site 2h) and P4mm (site 
lb). The values of ^ and 2 ;^ for the different models are listed 
in table AI. 


P4/mmm 


P4mm 


Fei 

(1/2,1/2,z) 

Fei 

(l/2,l/2,z) 

F 02 

(1/2,1/2,-z) 



Cui 

(l/2,l/2,z’) 

Cui 

(1/2,1/2,z’) 

CU2 

(1/2,1/2,-z’) 




A. Commensurate magnetic structure 

For the commensurate magnetic propagation vector 
kc = (1/2, 1/2, 1/2), kc is equivalent to -kc for both 
PA/mmm and PAmm. Moreover, the little group 
coincides with G for the two space groups. The gener- 
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TABLE IV. Labels of the generators of the space groups 
P4/mmm and P4mm, their symbols in Hermann-Maugin no¬ 
tation and their action on general positions within the unit 
cell. 


P4/mmm 

P4mm 

2z = 2 (0,0,z); (-X, -y, z) 

2z = 2 (0,0,z); (-X, -y, z) 

2y = 2 (0,y,z); (-x, y, -z) 

4+ = 4+ (0,0,z); (-y, x, z) 

4+ = 4+ (0,0,z); (-y, x, z) 

= m (x,0,z); (x, -y, z) 

i = -1 (0,0,0); (-X, -y, -z) 



ators of are listed in Tab. lEl and the characters of 
the irreps^ as calculated by the program BasiReps^ are 
displayed in Tab. [Vl 


TABLE V. Characters of the irreducible representations for 
the space groups P4/mmm and P4mm and the commensu¬ 
rate propagation vector kc = (1/2, 1/2, 1/2). Eor both space 
groups coincides with G. Only the characters correspond¬ 
ing to the generators of Gj^ are shown. 


P 4/171171171 

P4mm 

2z 2y 4+ i 

2z 4+ rrixz 



1 

1 

1 1 

rj 

1 

1 

1 

n 

1 

1 

1 -1 

ri 

1 

1 

-1 

n 

1 

1 

-1 1 

n 

1 

-1 

1 


1 

1 

-1 -1 


1 

-1 

-1 

n 

1 

-1 

1 1 

n 

-2 

0 

0 

n 

1 

-1 

1 -1 





n 

1 

-1 

-1 1 





n 

1 

-1 

-1 -1 





n 

-2 

0 

0 -2 






-2 

0 

0 2 






A reducible representation E can be constructed using 
the transformation properties of the magnetic moment 
components at the Fe sites. The decomposition of F in 
terms of irreps for Gj^ = PA/mmm is: 

r = r§ + r^ + rg + rjo (3) 

and for Gj^ = PAmm: 

r = r^2 + n ( 4 ) 


TABLE VI. Basis functions for the Fe atoms at the site (2h) of 
the space group P4/mmm for the commensurate propagation 
vector kc = (1/2, 1/2, 1/2). They also apply to the Cu atoms, 
which occupy the same crystallographic site with a different 
z coordinate. 


P4/i7imi7i (kc) 

IR 

Basis 

Fei 

Fe2 


vectors 

(1/2, 1/2, z) 

(1/2, 1/2, -z) 

r? 

- 

- 

- 

ri 


Re: [0,0,1] 

[0,0,-l] 


- 

- 

- 

r2 

- 

- 

- 

r§ 

- 

- 

- 

r§ 

- 

- 

- 

r? 


Re: [0,0,1] 

[0,0,1] 

ri 

- 

- 

- 

ri 

V® 

Re: [1,0,0] 

[-1,0,0] 



Re: [0,-l,0] 

[0,1,0] 

r?o 


Re: [0,1,0] 

[0,1,0] 


vr 

Re: [1,0,0] 

[1,0,0] 


The basis vectors obtained using the projector operator 
method for each one of the irreps appearing in the decom¬ 
position are listed in Tables [VT] and IVIIIl For PA/mmm 
the Fe magnetic moments can be along the c axis (F^ 
and Fy) or within the ab plane (Fg and Ffg)- For F^ and 
Fg the coupling between the moments at the Fei and Fe 2 
sites is antiferromagnetic (AFM) whereas for Tj and F^q 
it is ferromagnetic (FM). 

In the case of non-centrosymmetric PAmm the Fe mo¬ 
ments can only be along the c axis for F^ . For Fg , which 
is two-dimensional and complex, the orientation of the 
Fe moment can’t be obtained in a direct way. However, 
by imposing the condition that the magnetic moments 
should be real, it is easy to show that any linear com¬ 
bination of basis vectors gives rise to magnetic moments 
lying within the ab plane. To be noted is that for both 
irreps^ the phase between the moments at sites 1 and 2 
can not be derived from symmetry arguments and has to 
be obtained from the fit of the experimental data. 

Although symmetry arguments require the Gu/Fe 
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YBaCuFeOs, 230K DMC, = 2.458 A 




P4/nnnnnn P4nnnn 

ji] r/ + r,oC r/ + r,c 
jM]r/ + rgC r,c + r3c 


FIG. 10. (Color online) Left panel: Selected regions of the RT neutron powder diffraction pattern of YBaCuFeOs measured at 
DMC with A = 2.458A. Red crosses indicate experimental data, and continuous lines the fits carried out using the models for 
the magnetic structure schematized in the right panel. The irreps describing the different magnetic orders in the sapce groups 
P4/mmm and P4mmm are also indicated. 


TABLE VII. Characters of the irreducible representations for 
the space groups P4/mmm and P4mm and the ICM propaga¬ 
tion vector ki = (1/2, 1/2, fe). For both space groups = 
P4mm. Only the characters corresponding to the generators 
of Gj^ are shown. 


P 4/171171171 

P4i7im 

2z 4+ nixz 

2z 4+ mxz 


ri 

1 

1 

1 

ri 

1 

1 

1 


1 

1 

-1 

ri 

1 

1 

-1 

^ 3 

1 

-1 

1 

pZ 

1 

-1 

1 

ri 

1 

-1 

-1 

ri 

1 

-1 

-1 

r| 

-2 

0 

0 

ri 

-2 

0 

0 


magnetic moments to be either along the c axis or within 
the ab plane, no satisfactory fits could be achieved with 
these restrictions. Combining two irreps, two collinear 
arrangements and a commensurate helix with pitch = 
^ give reasonably good agreements. They are shown in 
Fig. [IQ] together with the associated irreps for each space 
group. As mentioned in section V, the three models are 
indistinguishable if Cu and Fe are sitting at z=l/4 and 
3/4. Using the actual atomic coordinates it is however 
possible to distinguish them. This is shown in Fig. [TOl 
where the results obtained for the three models are dis¬ 
played. For the Fe/Cu distribution we used the Model 
4 of Table i although identical conclusions were ob- 


TABLE VIII. Basis functions for Ee atoms at the site (2h) of 
the space group P4/mmm for the IGM propagation vector 
= (1/2, 1/2, kz). The table applies to each one of the two 
orbits in which the site (2h) splits for k^. The same table is 
obtained for the site (lb) of the space group P4mm with both 
ke = (1/2, 1/2, 1/2) and k, = (1/2, 1/2, kz) 


P4/i7imm (ki) 

P4i7im (kc and k^) 

IR 

Basis 

Fei 


vectors 

(1/2, 1/2, z) 

ri and Fi 

- 


F^ and ri 


Re: [0,0,1] 

Fs and ri 

- 


ri and Fi 

- 


ri and Fi 


Re: [1,0,0] 
Im: [0,-1,0] 


Vi 

Re: [-1,0,0] 
Im: [0,-1,0] 


tained for the other disordered arrangement (Model 2) 
and sligtly different values of the Fe and Cu magnetic 
moments. 
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Phase between sites 1 and 2 (degrees) 

371/5 271/5 



FIG. 11. (Color online) Upper pannels: variation of Chi2 as a function of the phase between the magnetic moments at the 1 
(blue) and 2 (red) sites. Lower panels: fits of the 1.5K data recorded on DMC with A = 2.45A using some representative values 
of the phases. Red crosses: experimental data. Black continuous lines: calculated patterns. The three rows of ticks indicate 
the positions of Bragg reflections of (up to down) the nuclear, the AFM commensurate and the AFM ICM phases. Particularly 
unstable fits (’’spikes” in a) marked as (2) and (5)) are obtained for values of the phase giving rise to spirals with uniform pitch 
(27r/5 and Ttt/S). 


B. Incommensurate magnetic structure 

For the ICM magnetic structure = (1/2, 1/2, kz) 
is not equivalent to — k^, neither for PA/mmm nor for 
PAmm. In the first case the possible magnetic structures 
are described by two arms of the star (k^ and —k^). For 
P4mm, whose star has only one arm (k^), we need to 
add the propagation vector k = — k^ for the description 
of magnetic structures with real magnetic moments. We 
note also that the little group = PAmm is the same 
for both PA/mmm and PAmm. 

The characters of the irreps of Gj^ calculated by the 
program BasiRepsi^ are displayed in Tab. IVIIl The de¬ 
composition of the reducible representation F in terms of 
irreps for Gj^ = PAmm is: 

r = q + ri. (5) 

The basis vectors associated with the irreps appearing 
in the decomposition are listed in Tab. Emi As in the 
case of the commensurate magnetic structure, the Fei 
moment can only be along the c axis for F^ or in the 
ab plane for F 5 . The phase between the moments of 
the ions sitting at the upper and lower pyramids can 
not be derived from symmetry arguments and has to be 
obtained from the fit of the experimental data. 

Fig. [TT] shows the evolution of the goodness of the fit 
as a function of the phase (j) between the magnetic mo¬ 
ments at the 1 (blue) and 2 (red) sites (note that the 


phase between two blue (two red) sites along c is given 
by the propagation vector). Only two values of 0 can re¬ 
produce the neutron powder diffraction data: (j) = tt and 
(j) = 7r/5. The first one involves the preservation of the 
AFM coupling between the red and blue sites without 
connecting oxygen observed in the commensurate phase 
and the loss of the FM coupling within the bipyramids. 

In the second, the moments within the bipyramids 
form an angle of 7r/2 as in the commensurate helix of 
Fig. [To] whereas the perpendicular arrangement is lost 
between the metallic sites without linking oxygen. In 
view of the smooth evolution between the commensurate 
and ICM structures, the first model (0 = tt) looks more 
plausible. Note also that it implies the existence of a less 
robust magnetic coupling within the bipyramidal units. 

We tried models based on variable-moment sinusoids 
and const ant-moment helices, either longitudinal, trans¬ 
verse or inclined with respect to the direction of ICM 
modulation vector. As for the commensurate magnetic 
structure, those with moments either along c or within 
the ab plane gave very poor agreements. The best re¬ 
sults were obtained for inclined circular helices, although 
equally good refinements could be obtained with a vari¬ 
able moment sinusoid and values of the Fe and Cu mo¬ 
ments slightly larger. In view of the observation of po¬ 
larization only by entering the ICM phase, the helical 
solutions were retained. 
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C. Calculation of exchange couplings 

Exchange coupling constants were extracted from ab 
initio calculation using collinear spin states (with spin- 
orbit coupling switched off). We applied the procedure 
described below^ assuming that all magnetic contribu¬ 
tions to the energy arise from the Heisenberg hamiltonian 
(Eq. 1 in the text). 

To extract the exchange coupling, Jab^ between ion A 
and ion B in the super cell, we calculate the energy for 
the following states: 1) all the spins in the unit cell are 
in the collinear commensurate spin order (i.e. the spin of 
ion A, Sao: is either parallel or antiparallel to the of ion 
B, Sbo, according to the commensurate spin order), 2) 
spin of ion A is flipped from its direction in the previous 
structure, 3) spin of ion B is flipped from its direction in 
configuration 1, 4) spin of ions A and B are both flipped. 
Assuming the magnetic ordering to be in the 2 : direction 
(we note that as spin-orbit coupling is not considered all 
directions are equivalent) and according to Eq. 1, the 
energies per supercell of these states are 

El =nJ abS%qS^q + ^aS^q + + e 

E2 = —bJabSao^bo ~ ^a^ao + e 

E3 = -uJabS^qS^o + ^aS^q - + e 

E/i =TiJ abSao^BO ~ ^A^ao ~ ^B^bo + ^5 

(6) 

where, hi is the effective field generated by all the spins 
in magnetic sublattices different from A or B and e is the 
sum of all the contributions (magnetic and non-magnetic) 
not involving spins at sublattices A and B. Moreover, n 
is the multiplicity of equivalent bonds connecting ion A 
to ion B when periodic boundary conditions are taken 


into account (e.g., using the four formula unit u.c. and 
the notation in Tab. I: n(Ji, 2 ) = 4, n{Ji^s) = 1 and 
^(-^ 3 , 5 ) = !)• Once the energies Ei are calculated ab- 
initio, Eqs. [6] become a system of equations in the un¬ 
known values Jab^ hA, hB and e which is solved to ex¬ 
tract Jab- 

Additionally, we note that, once further-neighbor cou¬ 
plings are considered and the unit cell used is not large 
enough, due to periodic boundary conditions, ions A and 
B might be connected by two inequivalent bonds. This 
happens, for example, in configuration a) when next- 
nearest-neighbor couplings along c are considered. In 
such a case, the NNN coupling Ji ^5 is not equivalent to 
NNN coupling Ji+c ,5 as the first exchange goes through a 
Cu-Cu-Ee pathway while the second goes through Cu-Ee- 
Ee pathway. It is not possible to extract these exchange 
couplings using the four formula units u.c. in Eig. as, 
for such supercell, the above method would yield only the 
average value of the two inequivalent couplings. There¬ 
fore, to extract separately these exchanges, it is necessary 
to double the supercell along the c direction. Indeed, to 
separate inequivalent NNN couplings for configuration a, 
c and d, a supercell doubled along the c direction to¬ 
gether with a T-centered k-point grid of size 8x8x2 was 
used in the calculations. 
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